Superconducting qubits can be coupled and addressed as trapped ions 



Yu-xi Liu, 1 - 2 L. F. Wei, 12 J. R. Johansson, 2 J. S. Tsai, 1 ' 2 ' 3 and Franco Nori 1 ' 2 ' 4 

'CREST, Japan Science and Technology Agency (JST), Kawaguchi, Saitama 332-0012, Japan 
2 Frontier Research System, The Institute of Physical and Chemical Research (RIKEN), Wako-shi, Saitama 351-0198, Japan 
3 NEC Fundamental Research Laboratories, Tsukuba, Ibaraki 305-8051, Japan 
4 Center for Theoretical Physics, Physics Department, Center for the Study of Complex Systems, 
The University of Michigan, Ann Arbor, Michigan 48109-1040, USA 
(Dated: February 2, 2008) 

Exploiting the intrinsic nonlinearity of superconducting Josephson junctions, we propose a scalable circuit 
with superconducting qubits (SCQs) which is very similar to the successful one now being used for trapped ions. 
The SCQs are coupled to the "vibrational" mode provided by a superconducting LC circuit or its equivalent 
(e.g., a SQUID). Both single-qubit rotations and qubit-LC-circuit couplings/decouplings can be controlled by 
the frequencies of the time-dependent magnetic fluxes. The circuit is scalable since the qubit-qubit interactions, 
mediated by the LC circuit, can be selectively performed, and the information transfer can be realized in a 
controllable way. 



PACS numbers: 03.67.Lx, 74.50.+r, 85.25.Cp 

I. INTRODUCTION 

Superconducting quantum circuits with Josephson junc- 
tions are currently studied for their potential applications in 
quantum information processing |1]. Quantum coherent os- 
cillations and conditional gate operations have been demon- 
strated using two coupled superconducting charge qubits [2, 
ID]. For a circuit with two coupled flux qubits, spectroscopic 
measurements show that it acts as a quantized four-level sys- 
tem J3l ■ Further, entanglement has been experimentally veri- 
fied in coupled flux ^ and phase Jfjl 0, H qubits. 

A major challenge for superconducting qubits (SCQs) is 
how to design an experimentally realizable circuit where the 
couplings for different qubits can be selectively switched on 
and off, and then scaled up to many qubits. Although two- 
qubit gates can be generated (see, e.g., Ref. Q9[|) with always- 
on interbit couplings, it is still verydifficult to scale up ex- 
perimental circuits iB La , [j, La B , B lal- Theoretically, the 
circuits (e.g., Refs. fifl MM M M S LH LH LUl) 
can be scaled up via a common data bus (DB). The DB 
modes are virtually excited (e .g., Refs . ifToL [Till ) or excited 
(e.g., Refs. d M El M [H \M). In the former 
case ifToL [Till , the effective qubit couplings can be switched 
on and off by changing the magnetic flux through the circuit 
within nanoseconds, which is a challenge for current experi- 
ments. In the later case lfT2l IT3I IT4I IT5I Il6l IT7ll , the qubit and 
the DB are coupled or decoupled when they have the same 
(resonant) or different (non-resonant) frequencies, realized via 
a sudden non-adiabatic change of either the qubit or the DB 
eigenfrequency. This introduces additional noise. 

The superconducting Josephson junction is a key building 
block of superconducting quantum circuits. Nonlinearity is its 
intrinsic characteristic. This nonlinearity can be used to adjust 
the inter-qubit couplings II 19U2(X fell 12211 by changing the cur- 
rent bias of the coupler, and thus cancelling the direct mutual 
inductance between the qubits. It can also be used to realize 
the switchable coupling between two inductively coupled su- 
perconducting flux qubits via a variable-frequency magnetic 
flux ll23ll . Also, recently, the level quantization of the LC cir- 



cuit has been experimentally demonstrated 11241,125m . 

Combining the variable-frequency-controlled coupling ap- 
proach i23ll and experimental achievements of the quantum 
LC circuit l24l l25ll . we now study a different approach to 
realize scalable SCQs via a common DB, which is either a 
quantum LC circuit or its equivalent, modeled by a harmonic 
oscillator 11241 12511 . The equivalent LC circuits can be either 
a cavity field (e.g., (l2l Hi) realized by, for instance, a 
one-dimensional transmission line resonator [16] or a super- 
conducting loop with Josephson junctions (e.g., a dc -based 
SQUID). More significantly, all SCQs can work at their opti- 
mal points when the data bus is a superconducting loop with 
Josephson junctions (this is not the case with standard LC 
DBs). In our this novel approach, the individual properties 
(e.g., eigenfrequencies) of the DB and SCQs are always fixed, 
but the SCQ-DB couplings can be conveniently controlled by 
changing the frequencies of the applied fame-dependent mag- 
netic fluxes (TDMFs). This is promising, because it is of- 
ten easier to produce fast and precise frequency shifts of the 
radio-frequency control signal in experiments, as opposed to 
changing the amplitude of the dc signal. 

We should point out that in our proposal the quantum LC 
circuit or its equivalent has to be excited when the information 
is transferred from one qubit to another; therefore it is an ac- 
tive element, not a passive one which is just virtually excited. 
Our proposal can be essentially reduced to the one used for 
trapped ions 112611 - The SCQs are coupled to the "vibrational" 
mode provided by a superconducting LC circuit or its equiv- 
alent (e.g., a SQUID). Both single-qubit rotations and qubit- 
LC-circuit couplings/decouplings can be controlled by the fre- 
quencies of the time-dependent magnetic fluxes. It means that 
SQCs can be coupled and separately addressed similarly to 
trapped ions. This similarity is significant because trapped 
ions ll26tl are further ahead, along the quantum computing 
roadmap. It is important to stress that our theoretical model 
can well explain the blue/red sideband excitations which have 
been ex peri mentally implemented in superconducting qubit 
circuits 11241 |27j| . 
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FIG. 1 : (Color online) The Zth flux qubit with three junctions is cou- 
pled to an LC circuit in (a) or a dc-biased SQUID with biased current 
lb in (b) by the mutual inductance {1 = 1, 2). (a) An externally 
applied magnetic flux through the Zth qubit loop includes a dc $1'' 
term and ac $1'' (t) term controlling the coupling in (a). The currents 
through the first qubit, second qubit, and LC circuit in (a) [or SQUID 
loop in (b)] are l">, j' 2 - 1 , and / respectively, (b) However, when a 
dc-biased SQUID forms an equivalent LC circuit, the SQUID-qubit 
couplings are controlled by a TDMF, &e(t) = Acos(ui c t), through 
the SQUID loop. The TDMF is added to the nonlinear qubit in (a) 
and to the nonlinear SQUID loop in (b). The configuration in (b) 
is significantly better, because both qubits can work at the optimal 
point / = 1/2. 



H. MODEL 

We study three-junction flux qubits (e.g., l28l l29ll ). As 
shown in Fig. Q] we consider the simplest circuit, where two 
flux qubits are coupled to a DB: either an LC circuit or a su- 
perconducting loop with junctions (e.g., a dc-biased SQUID). 
Without loss of generality and for simplicity, the DB is here 
assumed to be an LC circuit with an inductance L and a capac- 
itance C. The mutual inductance between the Zth qubit and the 
LC circuit is (1 = 1, 2). The applied magnetic flux 
through the Zth qubit loop in Fig.QJa) is assumed to include a 

s(') 



static (or dc) magnetic flux $0 , and also a TDMF 

§£'>(*) = 4 cos(4 l h), 



(i) 



which controls the qubit-DB couplings. Neglecting the mu- 
tual inductance between the two qubits, the Hamiltonian can 
be written as 



H 
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with the current / and magnetic flux <ft = IL through the LC 
circuit loop. Considering a three-junction qubit, the single- 
qubit Hamiltonian Hi in Eq. (O should be 
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after neglecting the qubit self-inductance and constant terms 
Each junction in the Zth qubit has a 
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capacitance C$, phase drop ipf\ and supercurrent jjp 



JqP sin ipf', with critical current 1$ . The loop current of the 



Zth qubit is 



where 
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with the convention Cjg = aiC^l = ajCj^ , and ai < 1. 
The LC circuit can be modeled by a harmonic oscillator de- 
scribed by the creation operator 
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and its conjugate a, with frequency lo = l/\/LC. Consider- 
ing the TDMF, the phase constraint condition [28] through the 
Zth qubit loop becomes 
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with the reduced bias flux 

/ = ($« - M«/)/$ . (8) 

Here, the bias / includes the flux M^'I, produced by the LC 
circuit. Thus, in the qubit basis, Eq. (|2]i becomes 
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after neglecting the constant terms. Here the Pauli operators 

of the Zth qubit are defined by 0+ — \ei){gi\, cr_ = |<7j)(ej|, 

and Oz = |ej)(ej| — \gi)(gi\- The computational basis states 

of the Zth qubit are defined 00, for $%\t) = 0, by the 
two lowest eigenstates, |0); = \gi), and = |e;), of Hi 
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with the two independent variables cp p 

= - <p®)/2. 
The first two terms in Eq. (O denote the free Hamiltoni- 
ans of both qubits and the LC circuit; u>q is the transition 
frequency of the Zth qubit. The always-on interaction Hamil- 
tonian between the Zth qubit and the DB in the third term of 
Eq. © is 



^ = (o t + o)(G I£ rW+h.c) 
with the coupling constant 
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The fourth term in Eq. (0 represents the interaction between 
the Ith qubit and its TDMF with the interaction strength 
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The fifth term of Eq. (0 is the controllable nonlinear inter- 
action among the Zth qubit, the DB, and the TDMF, with the 
coupling strength 
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This nonlinear interaction term between the Zth qubit, the DB, 
and the TDMF originates from the expansion of the loop cur- 
rent l"> of the Zth qubit in Eq. (0 to first order on the small 

reduced flux $i^(t)/$o v ' a me ph ase constrain condition in 
Eq. (0. Above, the TDMF $i (t) equals zero when calcu- 
lating the coupling strengths Gi, A;, and il;. That is, Iq and 
ljp are supercurrents through the loop and the third junction 
respectively when $e ' (t) = 0. 



III. SWITCHABLE INTERACTION BETWEEN QUBIT 
AND DATA BUS 

We find that the Hamiltonian (0 can be reduced to the one 
used in trapped ions ll26ll if the always-on interaction terms 
HQ can be neglected. This approximation is valid |l§] during 
the TDMF operations, in the large detuning regime between 
any qubit (e.g., Ith qubit) and the DB 
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which can be achieved when the circuit is initially fabricated. 
Thus, neglecting the always-on coupling H^ t between the 
data bus and the qubits, the Hamiltonian (0 is reduced to 
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which now has the same form as the one used for quantum 
computing with trapped ions, in the standard Lamb-Dicke 
limit (see, e.g., [26]). 

Therefore, the essential difference between our Hamilto- 
nian in Eq. (0 and the one used for experiments lfl6ll24ll25ll is 
that: (a) the nonlinear coupling between the data bus, qubits, 
and the classical field in Eq. (0 is very important for the su- 
perconducting case. Using this term, we can explain the side- 
band transitions in the experiments J24ll25ll : (b) the always-on 
coupling ifS between the qubits and the data bus should be 
negligibly small in our proposal. Our theoretical model is in 



contrast with those in Refs. 11241 12511 where: (a) there is no 
nonlinear coupling between the data bus, qubits, and the clas- 
sical field; (b) the always-on Hamiltonian H^ t could not be 
neglected. That is, in Refs. l24l l25ll . the Hamiltonian is just 
the usual Jaynes-Cummings model which cannot be directly 
used to explain the sideband excitations, especially for the ex- 
perimental results in Refs. 11241 12711 . 

Analogous to the case of trapped ions, in our proposed de- 
vices, three-types of dynamical evolutions {carrier process; 
red sideband excitation; and blue sideband excitation) can be 
produced by the TDMF using the frequency -matching (reso- 
nant) condition and neglecting all fast oscillating terms. These 
three dynamical evolutions can be described as follows. 



(i) If w, 



(i) 



, the qubit and the DB evolve indepen- 



dently in the large detuning condition. The external flux 
$e (t) is only used to separately address the Ith qubit rota- 
tions. These rotations are governed by the Hamiltonian 



h.C. 



(16) 



in the interaction picture and using the rotating-wave approxi- 
mation (RWA) (also for the and Hi 1 ' shown below). This 
is the so-called carrier process in the trapped ions approach. 



(ii) If the frequencies satisfy the condition lu 



,(0 



then the $e (i) assists the Ith qubit to couple resonantly with 
the DB. This is the red sideband excitation, governed by the 
Hamiltonian 



(17) 



(iii) In the blue sideband excitation, the frequencies satisfy 



the condition luX 1 ' = wi'' 



H[ l> = n ia a^> +h.c 



u>, with the Hamiltonian 



(18) 



Based on above discussions, it can be easily found that our 
derived Hamiltonian in Eq. (0, reduced to Eq. (15[ , can natu- 
rally explain experimental results on the sideband excitations. 
For example, in Ref. [27], the qubit and the DB frequencies 
are 14 GHz and 4.3 GHz, respectively; the frequency u> c for 
the red or blue sideband excitation is 9.7 GHz or 18.32 GHz. 
However, the Jaynes-Cummings model cannot be used to ex- 
plain these experiments. The qubit-DB couplings/decouplings 
can be controlledhy appropriately selecting the a4 of $e (i) 
to match/mismatch the above frequency conditions of the side- 
band excitations. 



IV. SINGLE- AND TWO-QUBIT GATES 



For the Ith qubit, the carrier process described by Hc^ can 
be used to perform the single-qubit rotation 
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Here, pi = | A; \t/H depends on the Rabi frequency | A; \/K and 
duration r; <j>i is related to the phase of the TDMF applied to 
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the Zth qubit. For example, the phases 0; = and 0/ = 37r/2 

correspond to the rotations (Pi) and R$ {Pi), about the x 
and y axis, respectively. Thus, any single-qubit operation can 
be realized by a series of B2 ) (Pi) and R { y l) (Pi) rotations with 
well-chosen different angles Pi. 

Two-qubit gates can be obtained using two qubits interact- 
ing sequentially with their DB as in Ref. fl26ll . There, the 
controlled phase-flip and the controlled-NOT gates can be ob- 
tained by three and five steps, respectively. Here, we only 
discuss the difference between our proposal and the one used 
for trapped ions. In our proposed circuit, the ratio |G;|/A; 
cannot be infinitely small. Then, the uncontrollable qubit- 
DB interaction £/S needs to be considered by the effective 
Hamiltonian tl30ri 
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when the kh qubit is not addressed by the TDMF. After in- 
cluding this effect, three pulses (successively applied to the 
first, second, and first qubits) with durations t\, t 2 , and t% 
(used to perform a controlled phase-flip gate in Ref. 12611 ) will 
result in a two-qubit gate U two - This can be expressed as 




Uu 



in the two-qubit basis {|<?i)|<? 2 ) 
with the parameters 

2|G 2 | 2 
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(22) 
(23) 
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The two-qubit gate [/two in Eq. (fJB is just a controlled phase- 
flip gate when the large detuning condition |G;|/A; ~ is 
satisfied. Moreover, any quantum operation can also be real- 
ized by combining the two-qubit gate E/two with other single- 
qubit operations. 



V. ENTANGLEMENT AND STATE TRANSFER 

We now consider two different external fields satisfying 
frequency-matching conditions, e.g., for the red sideband ex- 
citation, which are simultaneously applied to the two qubits in 
Fig. 03 Then, in the interaction picture and the RWA, the in- 
teraction Hamiltonian in Eq. (O, between the LC circuit and 
the two qubits, becomes 



= ^(OaV+Lc). 



(25) 



FIG. 2: (Color online) Plots of the /-dependent always-on qubit- 
DB coupling strength Gi (red curve) and the TDMF-controlled 
qubit-DB coupling strength (blue curve), rescaled by Ri — 

(27r/$ )M (,) E { jly/hu/2L. 



For simplicity, the coupling strengths between the LC circuit 
and different qubits are now assumed to be identical, e.g., 
fii = fi 2 = |il|e _ * e . If the LC circuit is initially prepared 
in the first excited state |1), then the wave-function of 
the whole system can be written as 



|*(t)) = -ie ie sm(V2Qt) [|ei)| 52 )|0) + | ffl )|e 2 }|0)] 

+ cos(V2m) |si>|fla)|l>. (26) 

When y/2Qt/H = tt/2, then the LC circuit is in the vacuum 
state |0) and a maximally entangled state between two qubits 
can be generated as 



I* 



!2 = -j=[\ei)\g 2 ) + |ffi)|e 2 )]. 



(27) 



i=i 



When adding one more qubit to Fig. [TJa) or (b), an un- 
known state = Pi\gi) + /3 2 |ei) in the first qubit can 
be transferred to the third one using the standard teleporta- 
tion procedure: i) a maximally entangled state |5 ,+ ) 2 3 = 
[|e 2 )|g3) + |g 2 )|e3)]/\/2 between the second and third qubits 
is prepared using the same method outlined above; ii) a CNOT 

(12) 

gate E/cnot is implemented for the first and second qubits 
(here, the second one is the target); iii) a Hadmard gate is im- 
plemented on the first one; iv) a simultaneous measurement, 
which can now be done experimentally [8], is performed on 
the first and the second qubits. The four different measured 
results {|ei,e 2 ), |ei,g 2 ), \gi,e 2 ), and |5i,.g 2 )} correspond to 
four outputs {l^i), |^2), |^3)> and \1jj3)} in the third qubit. 
The unknown state in the first qubit can be transferred to the 
third one when the measured result for the first and second 
qubits is \e±, e 2 ). However, appropriate gates (i.e., cri 3 \ 

(3) (3) 

and a z Vi ) need to be performed on the other three outputs 
mentioned above to transfer to the third qubit. 



5 



VI. EXPERIMENTALLY ACCESSIBLE PARAMETERS 

We now analyze the coupling constants related to the Ith 
qubit: (i) the always-on qubit-DB coupling strength 

G z oc( e/ |J %), (28) 
and (ii) the TDMF-controlled qubit-DB coupling strength 

nioc(ei|cos(2vJp + 27r/)|5i). (29) 

At the degeneracy point / = 1/2, the qubit potential is 
symmetric [29] and its ground and excited states have oppo- 
site parities; however, cos(2<y9 p + 2tt/) and the qubit loop cur- 
rent Iq 1 ^ have even and odd parities, respectively. Therefore 
fit = but Gi ^ when / = 1/2. Clearly, that fij = 
can be avoided by slightly shifting / away from the degener- 
acy point. The experiments on sideband excitations, e.g., in 
Refs. Il24ll27ll . were performed with f ^ 1/2. Moreover, the 
controlled phase-flip gate [26], requiring a transition from the 
ground state to the second excited state, also implies that the 
reduced bias Q flux / ^ 1 /2. 

Figure [2] shows the /-dependent coupling strengths G; and 
17;, rescaled by 

Ri = (2vr/$ )M (i) ^3 ) v /W 2i - (30) 

As in Ref. l24ll . the Josephson energies of the two bigger junc- 
tions of the hh qubit are E^ l { = E^ = 225 GHz and the ratio 
between the small and big junction areas is a; = 0.76. The ra- 
tio between the Josephson energy Ej\ and the charge energy 

Ec^ of the Zth qubit is about 30.8. Using the qubit parameters 
listed above and also taking the amplitude Ai — $0/30 of the 
TDMF applied to the ith qubit, Fig.|2]shows that G; and 17; are 
comparable when / is away from 1/2; e.g., Gi « 0.0579 Ri 
and 17; rs 0.0224 R t when / = 0.49. The strength 17; 
can also be larger than the strength Gi in the range, e.g., 
0-47 </< 0.477. 

If the capacitance and inductance of the LC circuit is 
taken Il24ll as 12 pF and 250 pH, then the frequency u> of 
the LC circuit is about 2.9 GHz. When the mutual induc- 
tance between the Ith qubit and the LC circuit is taken 
as 20 pH, then G; w 37.6 MHz and 17; w 14.6 MHz when 
/ = 0.49. The Ith qubit frequency computed is about 18 
GHz when / = 0.49. Therefore, the detuning between the 
Ith qubit and the LC circuit is A; w 15.1 GHz, and the ra- 
tio G1/A1 w 0.0015. Indeed, the always-on coupling G; is 
negligibly small when the Ith qubit works at / = 0.49 for 
measuring the sideband excitations. The phase corrections 9{ 
in Eq. (l2Tb should be very small with short operation times for 
those qubits when no $e (t) is applied. 

For the LC circuit, if its capacitance G and inductance L 
are assumed as ~ 1 pF and ~ 10 nH, respectively, then the 
LC circuit plasma frequency can be ~ 1.6 GHz. The linear 
dimension for the LC circuit can be < 1 cm. The estimated 
distance for a negligible mutual inductance between two near- 
est qubits is ~ 200 /im, and thus one DB can approximately 
interact with ~ 40 qubits. Of course, the larger L of the LC 



circuit could have a larger linear dimension (allowing, e.g., 
L ~ 100 nH), and then more qubits, here about 400, could 
interact with the LC circuit. In practice, the superpositions of 
the ground and excited states for an LC circuit decay on a time 
scale given by 1/ RC, here R is the residual resistance of the 
circuit and its radiation losses. 



VII. DISCUSSIONS AND CONCLUSIONS 

For flux qubits, the single-qubit states can be measured by 
using, e.g., either a tank circuit weakly coupled to the qubit [5] 
or a dc-SQUID [31]. If only a single-qubit measurement can 
be done at a time fl32l1 or simultaneous measurements can be 
done (e.g., as for phase qubits [33]), then any unknown quan- 
tum state can be reconstructed 11321 13311 and the information of 
the qubits can be read out. 

In our proposal, two crucial points are: (1) the qubit and 
the LC circuit data bus should initially have a large detuning, 
such that their always-on coupling is negligibly small when 
the TDMF-assisted qubit-DB coupling is implemented; (2) 
the nonlinearity of the JJs is essential to achieve our goal: 
the nonlinear coupling between these three: the qubit, DB, 
and TDMF. Based on these two requirements, the circuit can 
be modified according to different experimental setups, e.g., 
the LC circuit can be replaced by a superconducting loop 
with JJs (e.g., a dc -biased SQUID as in Fig. [lib)). Three- 
junction flux qubits can also be replaced by other qubits JT|], 
e.g., one- or four-junction flux qubits, phase qubits, or charge- 
flux qubits. Although the self-inductances of the qubits are 
neglected here, our method is still valid for the qubits with 
nonzero self-inductances [13411 . 

Our numerical calculations show that the TDMF-controlled 
coupling strength 17; is not large enough to realize very fast 
two-qubit operations when the DB is a simple LC circuit. 
In principle, this problem could be solved by using a super- 
conducting loop with Josephson junctions (e.g., dc-SQUID in 
Ref. IU9I1 ) as a data bus instead of a simple LC circuit. Thus, 
the TDMF can be applied to the DB loop and the qubit can 
work at the optimal point; the DB-qubit always-on coupling 
can be minimized to zero; the TDMF-controlled coupling- 
strength can be large enough to realize fast two-qubit oper- 
ations. A more detailed study on this issue will be presented 
elsewhere. 

In conclusions, using the nonlinearity of the supercon- 
ducting JJs, we theoretically explain the sideband excita- 
tions for qubits coupled to an LC circuit and show how to 
scale these to many qubits. In contrast to previous propos- 
als (e.g., lfl^[Ill[m[l7ll ). the properties (e.g., eigenfrequen- 
cies) of the qubits and the DB are fixed when processing either 
the resonant coupling or the non-resonant decoupling. Also, 
the qubit-DB couplings/decouplings are controlled neither by 
changing the magnetic flux through the loop nor by changing 
the eigenfrequencies of the qubits (or the DB). They are only 
controlled via the frequency shifts of TDMFs, which is much 
easier to achieve experimentally. 
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